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ABSTRACT 

We use new central stellar velocity dispersions and nuclear X-ray and Ha luminosities for the Palomar survey 
of nearby galaxies to investigate the distribution of nuclear bolometric luminosities and Eddington ratios for 
their central black holes (BHs). This information helps to constrain the nature of their accretion flows and the 
physical drivers that control the spectral diversity of nearby active galactic nuclei. The characteristic values of 
the bolometric luminosities and Eddington ratios, which span over 7-8 orders of magnitude, from Lbd <> 10 37 to 



1 , vary systematically with nuclear spectral classification, increasing 
-> transition objects — > low-ionization nuclear emission-line regions 



3 x 10 44 erg s" 1 and W^Edd ~ 10" 9 to 10 
along the sequence absorption-line nuclei 
— > Seyferts. The Eddington ratio also increases from early-type to late-type galaxies. We show that the very 
modest accretion rates inferred from the nuclear luminosities can be readily supplied through local mass loss from 
evolved stars and Bondi accretion of hot gas, without appealing to additional fueling mechanisms such as angular 
momentum transport on larger scales. Indeed, we argue that the fuel reservoir generated by local processes should 
produce far more active nuclei than is actually observed. This generic luminosity-deficit problem suggests that 
the radiative efficiency in these systems is much less than the canonical value of 0. 1 for traditional optically thick, 
geometrically thin accretion disks. The observed values of Lboi/^Edd, all substantially below unity, further support 
the hypothesis that massive BHs in most nearby galaxies reside in a low or quiescent state, sustained by accretion 
through a radiatively inefficient mode. 

Subject headings: black hole physics — galaxies: active — galaxies: nuclei — galaxies: Seyfert 



1. INTRODUCTION 

Simple considerations of the quasar population predict that 
massive black holes (BHs) should be common in a sizable frac- 
tion of present-day galaxies. From the integrated luminosity 
density of quasars, one can estimate that a typical L* galaxy 
should, on average, contain a waste mass of ~ 10 7 M Q locked 
up in a BH (Soltan 1982; Chokshi & Turner 1992). But why 
are the quasar remnants so quiescent? Active galactic nuclei 
(AGNs) with quasarlike luminosities are absent at z ~ pre- 
sumably because of the diminished fuel supply currently avail- 
able. For a canonical radiative efficiency of r\ = 0.1 appropri- 
ate for geometrically thin, optically thick accretion disks (see 
Frank et al. 1992), the BH has to consume 1-100 M Q yr" 1 in 
order to generate luminosities of Lboi = 10 n -10 13 L©, as typ- 
ically seen in quasars. This level of gas supply is difficult to 
sustain in nearby galaxies. On the other hand, accretion rates 
of M = 0.001-0.1 Mq yr" 1 do not appear implausible. Even 
if angular momentum transport on nuclear scales is ineffective 
in disk galaxies, this level of fueling can be supplied simply 
through local stellar mass loss (Ho et al. 1997c). Hence, for rj 
= 0.1, there ought to be many nuclei shining as AGNs with Lboi 
= 10 9 -10 n L Q . This is not observed. Only ~l%-5% of galax- 
ies contain bright Seyfert nuclei (e.g., Huchra & Burg 1992; 
Greene & Ho 2007a). In the case of giant elliptical galaxies ex- 
periencing cooling flows, we may expect even larger values of 
M and thus correspondingly larger luminosities, again contrary 
to observations (Fabian & Canizares 1988). 

The dilemma posed by the luminosity deficit in the nu- 
clei of nearby elliptical galaxies can be resolved by discarding 
the premise that massive BHs are ubiquitous in these systems 
(Fabian & Canizares 1988). But this proposition is no longer 
tenable in light of our current knowledge on the demography of 



central BHs based on direct dynamical searches (Magorrian et 
al. 1998; Ho 1999a; Kormendy 2004). Massive BHs appear to 
be a generic component of galactic structure in most, if not all, 
systems with a bulge. Consistent with this picture, low-level 
nuclear activity qualitatively resembling that of more luminous 
AGNs is found to be equally pervasive in nearby galaxies (Ho 
etal. 1997b; Ho 2008). 

Radiatively inefficient accretion flows (RIAFs; see Narayan 
et al. 1998; Quataert 2001 for reviews) provide an attractive 
framework for solving the luminosity-deficit problem. In the 
regime when the mass accretion rate onto the central BH is 
very low, the low-density, optically thin accreting medium can- 
not cool efficiently, and the accretion flow consequently puffs 
up into a quasi-spherical structure. Most relevant to the present 
discussion, RIAFs attain radiative efficiencies much below the 
canonical value of 0.1. RIAFs are characteristically dim. Opti- 
cally thin RIAFs are predicted to exist for accretion rates below 
a critical threshold of M cr ; t ps a 2 Medd ~ 0. lMEdd (Narayan et al. 
1998), where the Eddington accretion rate is defined by L E dd = 
T]M E ddC 2 , with 77 = 0.1 and L Edd = 1.26 x 10 38 (M BH /M Q ) erg 
s" 1 ; the Shakura & Sunyaev (1973) viscosity parameter is taken 
to be a fa 0.3 (Narayan et al. 1998). Fabian & Rees (1995) 
and Mahadevan (1997) invoke RIAFs to explain the dimness of 
elliptical galaxy nuclei. 

This paper discusses the luminosity-deficit problem for a 
large, well-defined sample of galaxies spanning a wide range 
of morphological types and representing all the major nuclear 
spectroscopic classes. We use X-ray and line luminosity mea- 
surements to constrain the accretion luminosities of the nuclei. 
A newly published catalog of central stellar velocity dispersions 
provides estimates of the BH mass through the Mbh-c* relation 
(Gebhardt et al. 2000; Ferrarese & Merritt 2000). The feeble 
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TABLE 1 



Bolometric Corrections for Low-Luminosity AGNs 



Galaxy 


Class 


D L 


Lboi 








C'x 


Reference 






(Mpc) 


(erg s" 1 ) 


Narrow 


Broad 


Total 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NOC 1097 


1,2/1,1 9 


14.5 


1.8 xlO 42 


529 


474 


250 


34.4 


1 


NGC 3031 


S1.5/L1.5 


3.6 


2.1 xlO 41 


300 


140 


95 


3.5 


2 


NGC 4203 


LI. 9 


9.7 


9.5 xlO 40 


431 


243 


156 


3.9 


3 


NGC 4261 


L2 


30.0 


1.7xl0 42 


355 




355 


14.2 


3 


NGC 4374 


L2 


16.8 


8.2xl0 41 


410 




410 


8.4 


3 


NGC 4450 


LI. 9 


16.8 


3.4xl0 41 


1064 


1136 


550 


8.9 


3 


NGC 4486 


L2 


16.8 


2.3 xlO 42 


389 




389 


17.7 


3 


NGC 4579 


S1.9/L1.9 


16.8 


9.9 xlO 41 


353 


380 


184 


7.7 


3 


NGC 4594 


L2 


9.2 


2.7xl0 41 


208 




208 


7.5 


3 


NGC 6251 


S2 


92.0 


8.0 xlO 42 


36 




36 


2.9 


4 


Arp 102B 


LI. 2 


96.6 


2.9 xlO 43 


546 


62 


55 


1.8 


5 


Pictor A 


L1.5 


140.1 


1.3 xlO 44 


5000 


241 


228 


2.7 


6 








Mean 


802 


382 


243 


9.5 










Standard Deviation 


1345 


360 


156 


9.2 










Median 


410 


243 


228 


7.7 





NOTE. — Col. (1) Galaxy name. Col. (2) Spectroscopic classification of nucleus. Col. (3) Adopted distance. Col. (4) 
Bolometric luminosity. Col. (5) Ratio of bolometric to narrow Ha luminosity. Col. (6) Ratio of bolometric to broad Ha 
luminosity. Col. (7) Ratio of bolometric to total (narrow + broad) Ha luminosity. Col. (8) Ratio of bolometric to X-ray 
luminosity in the 2-10 keV band. Col. (9) Reference for the Ha data: (1) Storchi-Bergmann et al. (1993), narrow-line 
luminosity corrected for an extinction of Ay = 2.46 mag; (2) Ho et al. (1996); (3) Ho et al. (1997a); (4) Shuder & 
Osterbrock (1981); (5) Halpern et al. (1996); (6) Carswell et al. (1984). Data for cols. (2)-(4) and (8) come from Ho 
(1999b) and Ho et al. (2000). 



nuclear activity in nearby galaxies can be comfortably sustained 
through stellar mass loss and spherical accretion of hot gas in 
the inner regions of bulges. We show that the inferred accretion 
rates lie well within the critical threshold of RIAFs. Finally, 
we discuss the physical connection between accretion rates and 
nuclear spectral types. 

2. OBSERVATIONAL MATERIAL 

Our analysis explicitly makes the following assumptions: (1) 
all galaxy bulges contain central BHs whose masses can be 
well estimated using the recently established correlation be- 
tween BH mass and bulge stellar velocity dispersion; (2) the 
present level of activity of the BHs manifests itself as AGN- 
like emission-line nuclei; and (3) the accretion luminosity can 
be extrapolated through the observed nuclear X-ray or optical 
line luminosities. 

We focus on the sample of nuclei in the Palomar survey of 
nearby galaxies, a magnitude-limited spectroscopic study of a 
nearly complete sample of 486 bright (B T < 12.5 mag), north- 
ern (5 >0°) galaxies (see Ho et al. 1997a, 1997b, and refer- 
ences therein). All of the galaxies have been assigned a nu- 
clear spectroscopic classification using a set of uniform crite- 
ria, into the following classes (see Ho et al. 1997a for details): 
absorption-line nuclei (A), H II nuclei (H), Seyfert nuclei (S), 



low-ionization nuclear emission-line regions (LINERs; L), and 
transition objects (T; LINER/H II composites). This paper con- 
siders all classes except the H II nuclei. The Palomar survey 
contains 277 galaxies classified as absorption-line (66), Seyfert 
(52), LINER (94), and transition nuclei (65). 

2.1. Nuclear Luminosities 

Our study requires bolometric luminosities for a well-defined 
sample of AGNs covering a wide range of power. Because 
AGNs emit a very broad spectrum, their bolometric luminosi- 
ties ideally should be measured directly from their broadband 
spectral energy distributions (SEDs). In practice, however, 
complete SEDs are not readily available for most AGNs, least 
of all for low-luminosity sources such as LINERs, which are 
most prevalent in nearby galaxies. The largest existing compila- 
tions of broadband SEDs for low-luminosity AGNs (Ho 1999b; 
Ho et al. 2000; Maoz 2007) contain only a limited number of 
objects. 

For this study, we circumvent this difficulty by using two 
measures of the nuclear power to estimate the AGN bolomet- 
ric luminosity, one based on the Ha emission line and another 
based on X-rays. Although the Ha luminosity comprises only 
a small percentage of the total power, its fractional contribu- 
tion to the bolometric luminosity, as shown below, turns out 
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to be fairly well defined. Moreover, unlike most other spec- 
tral windows, Ha measurements are readily available for large, 
relatively complete samples of nuclei. The Palomar database 
gives Ha luminosities measured through an aperture of 2" x4" 
centered on the nucleus, which corresponds to a linear scale 
of ^200 pc x 400 pc for a typical distance of 20 Mpc. As 
explained in Ho et al. (2003a), some of the Ha luminosities 
published in Ho et al. (1997a) have been updated with more 
accurate values from the literature. We will also make use of 
upper limits for the Ha luminosity of the absorption-line nuclei 
that were not given explicitly in Ho et al. (1997a); the limits 
were calculated from the equivalent-width detection limit of the 
survey in conjunction with the 6600 A continuum flux density 
measurements, assuming that the line has a typical full width 
at half-maximum (FWHM) of 250 km s" 1 . A supplementary 
list of Ha luminosities, including the upper limits, is given in 
Ho et al. (2003a). In total, Ha luminosities, or upper limits 
thereof, are available for 246 objects, which account for 80%, 
98%, 89%, and 91% of the nuclear classes A, S, L, and T, re- 
spectively. Thus, Ha measurements are available for the vast 
majority of the objects considered in this study. 

The bolometric correction for Ha, Ch q = L\, i/Lu a , can be 
obtained in one of two ways. For luminous, type 1 AGNs, it 
is often expedient to estimate Lboi from the optical continuum 
luminosity, frequently chosen at 5100 A in the recent litera- 
ture: Lboi = C5iooALa(5100). Now, the Ha luminosity corre- 
lates strongly with the optical continuum, including ALa(5100) 
(Greene & Ho 2005); the correlation is slightly nonlinear. 
Choosing C5100 = 9.8 (McLure & Dunlop 2004), Greene & Ho 
(2007a) obtain L bol = 2.34 x 10 44 (L Hce /10 42 erg s" 1 ) 86 erg s _1 . 
The conversion pertains to the entire Ha line, which in lumi- 
nous type 1 AGNs is dominated by the broad component. Be- 
cause of the reliance on the L Ka -\L\(5 100) correlation and the 
assumption that C5100 is constant, the relation between Lboi and 
Lu a is formally slightly nonlinear. It is unclear how robust this 
result is and whether it can be extrapolated toward lower lumi- 
nosities. For the luminosity range of interest to us, it introduces 
an uncertainty of a factor of ~ 2 into the bolometric correction. 
For L Ha = 10 38 , 10 39 , and 10 40 erg s" 1 , C Ha « 850, 615, and 
446, respectively. 

Alternatively, we can attempt to estimate Ch q empirically 
from the observed SEDs of low-luminosity AGNs, limited 
though they may be. We use the sample of 12 low-luminosity 
Seyferts and LINERs with broadband SEDs studied by Ho 
(1999b) and Ho et al. (2000; see also Ho 2002b) as a guide. The 
data, summarized in Table 1, show that the median value of Ch q 
ranges from 228 to 410, depending on whether we include only 
the narrow component of the line, only the broad component, 
or both. A reasonable compromise might be Ch q ~ 300 ± 100. 
Preliminary analysis of a more extensive sample of SEDs (L. C. 
Ho, in preparation) shows that the bolometric correction for the 
Hf3 line has a significant scatter, especially for low-luminosity 
sources. For sources with Eddington ratios below 0.1, and as- 
suming that on average Ha/H/3= 3.5 (Greene & Ho 2005), the 
median value of Ch q ~ 220 with an interquartile range of 160. 
This agrees reasonably well with the range of values obtained 
from our calibration sample in Table 1 . For concreteness, we 
will simply adopt Cho = 300; none of the main conclusions in 
this study depends critically on the exact value of the bolomet- 
ric correction. 

The Ha luminosities of the narrow-line objects are subject 
to a potential source of complication. From consideration of 



the relative strength of X-ray (2-10 keV) and Ha emission in 
different classes of nearby AGNs, Ho (2008) shows that type 2 
sources (Seyfert 2s, LINER 2s, and essentially all of the tran- 
sition objects) have a general tendency to emit excess optical 
line emission compared to their type 1 counterparts. He at- 
tributes this excess emission to extranuclear processes unasso- 
ciated with the active nucleus. If this interpretation is correct, 
then only a fraction of the narrow Ha emission should be in- 
cluded in the budget for the nuclear luminosity. Ho (2008) 
finds that the median ratio of Lx/Lh q is 7.3 for Seyfert Is and 
4.6 for LINER Is, whereas Seyfert 2s, LINER 2s, and transi- 
tion objects have significantly lower values of 0.75, 1.6, and 
0.41, respectively. Assuming, for concreteness, that the in- 
trinsic Lx/LH a ratio for the type 2 sources is equal to that of 
LINER Is, the Ha luminosities for Seyfert 2s, LINER 2s, and 
transition objects need to be reduced by a factor of 6.1, 2.9, and 
11.2, respectively. 

In light of the above source of ambiguity with the Ha emis- 
sion, it would be desirable to have an alternative handle on 
the nuclear luminosity. By far the most secure measure of nu- 
clear luminosity in AGNs comes from X-ray observations. Be- 
cause of the faintness of low-luminosity AGNs, the data must 
be of high enough angular resolution so that the nucleus can be 
cleanly separated from the host galaxy (e.g., Ho et al. 2001; 
Flohic et al. 2006). As reviewed in Ho (2008), a substantial 
fraction of the galaxies in the Palomar survey have now been 
observed in the X-rays. The Appendix gives a compilation 
of all pertinent X-ray measurements taken from the literature, 
as well as from new analysis of data taken from the Chandra 
public archives. The X-ray data are not nearly as complete as 
Ha. Nevertheless, X-ray luminosities, or upper limits thereof, 
are now available for 175 out of the 277 objects in the parent 
sample (63%), which account for 47% of the absorption nu- 
clei, 68% of the LINERs, 83% of the Seyferts, and 57% of the 
transition objects. The incompleteness and heterogeneous na- 
ture of the X-ray measurements make it difficult to rigorously 
assess selection effects. However, if observational biases exist, 
they should be in the direction of missing very faint sources, an 
effect that strengthens our main conclusions. 

As with the Ha data, determining the appropriate bolometric 
correction for the X-ray band is not trivial. Since the SEDs of 
AGNs vary strongly with accretion rate (Ho 1999b), we must 
abandon the usual practice of adopting a single correction factor 
based on the average SED of luminous quasars. Using, again, 
the small sample of low-luminosity AGNs with reliable broad- 
band SEDs, Table 1 shows that the median value of the bolo- 
metric correction in the 2-10 keV band is Cx = Lb i/L x ~ 8, 
where L x is the luminosity in the 2-10 keV band, corrected, to 
the extent possible, for absorption. The more extensive data set 
of L. C. Ho (in preparation) suggests a value larger by about 
a factor of 2: sources with Lboi/^Edd <^ 0.1 have a median 
Cx = 15.8, with an interquartile range of 9.6. Because low- 
luminosity AGNs tend to be "X-ray-loud" (Ho 1999b) their val- 
ues of Cx are significantly smaller than conventionally assumed 
for luminous sources (C x ~ 35; Elvis et al. 1994). This is con- 
sistent with the analysis of Vasudevan & Fabian (2007). For the 
present purposes, we will adopt Cx = 15.8, noting, as before, 
that factors of a few variation in the bolometric correction do 
not affect the main conclusions of this study. 

The uncertainties on Lboi are difficult to estimate. As dis- 
cussed in Ho et al. (1997a), the Ha fluxes in the Palomar survey 
have typical errors of 30%-50%, reaching 100% in the worst 
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FIG. 1 . — Distribution of (a) bolometric luminosity, Lb i, and (b) ratio of bolometric luminosity to the Eddington luminosity Lbol/^Eddi Ior all objects, Seyferts 
(L), LINERs (L), transition nuclei (T), and absorption-line nuclei (A). The hatched and open histograms denote detections and upper limits, respectively. The 
original sample is shown in blue, and the subsample restricted to Sab-Sbc (T = 2—4) Hubble types is shown in red. The bolometric luminosity is based on the 2—10 
keV X-ray luminosity, assuming i^ol = 15.8Z.x- 
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- Distribution of (a) bolometric luminosity, Lboh and (b) ratio of bolometric luminosity to the Eddington luminosity, Lbol/Z-Edd, for galaxies binned 
type. The hatched and open histograms denote detections and upper limits, respectively. The bolometric luminosity is based on the 2-10 keV X-ray 
assuming Lb i = 15.8Lx- 



cases. The largest source of uncertainty for the Ha-based lumi- 
nosities, however, comes from our still-tentative knowledge of 
Cho (factor ~ 2) and the amount of extranuclear contamination 
of the narrow-line emission (factor ~ 2 — 5, depending on the 
spectral class). In the X-ray band, we can be more confident 
that the flux is largely nuclear, and large-amplitude variability 



seems to be rather uncommon for the systems in question (Ho 
2008). Not all of the X-ray detections have sufficient counts for 
rigorous spectral fitting, but fortunately low-luminosity AGNs 
generally have small absorbing columns (Ho 2008). Still, at the 
moment we do not know Cx to better than a factor of ~ 2. We 
conservatively guess that the estimates of Lbol based on Ha 
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FIG. 3. — Distribution of (a) bolometric luminosity, Lboi, and (b) ratio of bolometric luminosity to the Eddington luminosity, Lbol/^Eddi f° r all objects, Seyferts 
(L), LINERs (L), transition nuclei (T), and absorption-line nuclei (A). The bolometric luminosity is based on the extinction corrected total (narrow + broad) Ha 
luminosity, assuming = 300LH a ■ The blue histograms show the distributions after correcting the narrow-line sources for extranuclear contamination (see text 
for details). The hatched and open histograms denote detections and upper limits, respectively. 



and X-rays measurements have uncertainties of 0.7 and 0.3 dex, 
respectively. 

2.2. Black Hole Masses 

The majority of the objects in our sample do not have di- 
rect, dynamically determined BH masses. We will estimate BH 
masses using the tight correlation between BH mass and bulge 
stellar velocity dispersion (the Mbh-c* relation: Gebhardt et al. 
2000; Ferrarese & Merritt 2000), as determined by Tremaine et 
al. (2002): 

log f — ) = (4.02 ±0.32) log f — — r ) +(8.13 ±0.06). 

B \M J K g V200kms-7 

(1) 

The intrinsic scatter of the above fit is estimated to be <;0.3 
dex. The fit uses a e , the luminosity-weighted velocity disper- 
sion measured within the effective radius of the bulge. Since 
we do not have measurements of a e for most of our galaxies, 
we use instead ctq, the central velocity dispersion. Gebhardt et 
al. (2000) have shown that in general oo ~ &e within a scatter 
of -10%. 

Ho et al. (2009) recently published a comprehensive, uni- 
form catalog of central velocity dispersions for nearly all of the 
galaxies in the Palomar survey. New stellar velocity dispersions 
were obtained for a total of 428 galaxies, and estimates for an- 
other 34 were obtained indirectly from the line width of [N II] 
A6583 using the calibration of Ho (2009). The typical uncer- 
tainties in the velocity dispersions range from —5% to 15%. 
An error of 10% in o. t introduces an uncertainty of —0.15 dex 
in log Mbh- We assume that Mbh has an uncertainty dominated 
by the scatter of the Mbh-c* relation, —0.3 dex. 



3. BOLOMETRIC LUMINOSITIES AND EDDINGTON RATIOS 

Figure 1 shows the distributions of bolometric luminosities 
and their values normalized with respect to the Eddington lumi- 
nosity, Lboi/^Edd- To avoid the possible complication of extranu- 
clear contamination in the Ha emission, we base the bolometric 
luminosities on the hard X-ray measurements. The statistics of 
the distributions are listed in Table 2. The four classes of nuclei 
comprise a sequence of increasing luminosity: A — > T — > L — > 
S. Whereas LINER and transition nuclei have very similar Ha 
luminosities (Ho et al. 2003a) — an effect that can be attributed 
to the Ha emission in transition objects being boosted by non- 
nuclear sources (Ho 2008) — there is no doubt that in the hard 
X-ray band LINERs are more luminous than transition objects 
(median L bo i = 3.0 x 10 40 vs. 6.5 x 10 39 erg s" 1 ). Both LINERs 
and transition objects, in turn, are less powerful than Seyferts 
(median L\, \ = 2.2 x 10 41 erg s -1 ). This systematic trend be- 
comes even more sharply defined when we consider the Ed- 
dington ratios. As with Lboi, the median value of Lboi/^-Edd for 
LINERs (6.0 x 10~ 6 ) is a factor of 4 larger than for transition ob- 
jects (1.5 x 10~ 6 ), both being 1-2 orders of magnitude smaller 
than in Seyferts (1.1 x 10~ 4 ). Notably, all galactic nuclei in the 
Palomar sample are highly sub-Eddington systems. 

The distribution of Lboi is broadly similar for galaxies of dif- 
ferent morphological types (Fig. 2a). By contrast, Lboi/^Edd in- 
creases mildly, but systematically, from early-type to late-type 
galaxies (Fig. 2b). Since the various classes of emission-line 
nuclei in the Palomar survey are hosted by slightly different 
Hubble types (Ho et al. 2003a), it would be of interest to ex- 
amine the trends in Lboi and Lboi/Ledd after factoring out the 
dependence on Hubble type. This is illustrated by the red his- 
tograms in Figure 1, where we now restrict the comparison to 
galaxies with morphological types Sab-Sbc (T = 2-4), a subset 
that, as shown in Ho et al. (2003a), has statistically identical 
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TABLE 2 



Sample Statistics 



Sample 


N 






Lboi (erg s 1 ) 








ibol/iEdd 






Mean 


Error 


Median 


Mean 




Error 




Median 




All (all types) 


175 


37 


8.5xl0 42 


4.2 xlO 42 


2.0 xlO 40 


2.3x10" 


-3 


1.2x10" 


-3 


4.0x10" 


-6 


S (all types) 


43 


3 


3.4xl0 43 


1.7xl0 43 


2.2xl0 41 


9.2x10" 


-3 


4.9x10" 


-3 


1.1x10" 


■1 


(Sab-Sbc) 


15 





5.9xl0 43 


4.3 xlO 43 


7.2 xlO 41 


2.3x10" 


-2 


1.3x10" 


-2 


2.5x10" 


-4 


L (all types) 


64 


12 


6.1xl0 41 


2.3xl0 41 


3.0 xlO 40 


3.0x10" 




1.0x10" 




6.0x10" 


-6 


(Sab-Sbc) 


16 


2 


1.9xl0 41 


8.9 xlO 40 


8.3 xlO 40 


3.3x10" 




1.2x10" 


-5 


6.9x10" 


-6 


T (all types) 


37 


10 


3.0xl0 4 " 


9.7 xlO 39 


6.5 xlO 39 


5.1x10" 




3.7x10" 


-5 


1.5x10" 


-6 


(Sab-Sbc) 


17 


6 


3.5 xlO 40 


1.8 xlO 40 


3.4 xlO 39 


3.1x10" 




2.5x10" 


-5 


1.9x10" 


-6 


A (all types) 


31 


10 


1.3xl0 41 


5.8 xlO 40 


1.8 xlO 39 


3.1x10" 


-6 


1.0x10" 


-6 


2.2x10" 


-7 


E 


41 


10 


4.6xl0 41 


2.0xl0 41 


1.7X10 40 


1.2x10" 


-5 


5.8x10" 


-6 


1.2x10" 


-6 


SO-SO/a 


51 


12 


8.1xl0 42 


5.3 xlO 42 


1.2 xlO 40 


7.1x10" 


-4 


4.7x10" 


-4 


4.7x10" 


-6 


Sa-Sbc 


64 


9 


1.5xl0 4;i 


1.1 xlO 43 


3.1 xlO 40 


5.5x10" 


-3 


3.3x10" 


-3 


9.2x10" 


-6 


Sc-later 


19 


6 


6.1xl0 42 


5.7xl0 42 


3.6 xlO 39 


5.9x10" 


-4 


3.7x10" 


-4 


4.2x10" 


-6 



a The bolometric luminosities are based on the 2-10 keV X-ray luminosity, assuming Lboi = 15.8Lx. Statistics for 
subsamples containing upper limits (whose number is denoted by N u ) computed using the Kaplan-Meier product- 
limit estimator (Feigelson & Nelson 1985). Our sample contains four galaxies (NGC 1275, 4261, 4374, and 4486) with 
strong FR I-type (Fanaroff & Riley 1974) radio jets, as noted in Table 3. The nuclear X-ray luminosities of these 
objects might be partly contaminated by emission from the jet component. Excluding these four galaxies does not 
significantly change the statistics of this table. 



distributions of morphological types for all three AGN sub- 
classes. The overall trends of the parent sample are preserved 
(see also Table 2). 

Although in this study we give preference to the X-ray lu- 
minosities over the Ha luminosities because of concerns over 
extranuclear line contamination, the Ha data have the advan- 
tage of being uniform and nearly complete. Figure 3 repeats the 
analysis of Figure 1, but now using bolometric luminosities de- 
rived from Ha. The black histograms show the Ha luminosities 
as observed, converted to Lboi assuming a bolometric correction 
of Cua = 300. The blue histograms plot the same data with a 
statistical correction for extranuclear line emission applied to 
the Seyfert 2s, LINER 2s, and transition objects based on their 
observed L x /L Hq ratio (see § 2. 1). Not surprisingly, the overall 
trends seen in Figure 1 are well mirrored in Figure 3, but now 
they are delineated with better statistics. 

Finally, we turn to the distribution of Mbh and Lboi vs. 
Lboi/Ledd (Fig. 4). In these diagrams, we have marked the var- 
ious subclasses of nuclei, and we have included the large sam- 
ple of z < 0.35 high-luminosity AGNs (Seyfert 1 nuclei and 
low-redshift quasars) selected by Greene & Ho (2007a) from 
the Sloan Digital Sky Survey (SDSS). Greene & Ho derived 
BH masses and bolometric luminosities using the width and 
strength of the broad Ha line. We point out several salient fea- 
tures. 

1 . Considered collectively, there is no dependence of M SH 
on Lboi/LEdd^ at a given value of Mbh, which mostly fall 
in the range 10 6 — 10 9 M©, Lboi/LEdd spans ~ 6 orders 
of magnitude within the Palomar sample and <~ 8 orders 
of magnitude if the sample of luminous sources is in- 
cluded. 

2. Within each class of emission-line objects, especially 
among LINERs and transition nuclei, there is a loose 



inverse correlation between M B h and L\, \/L E dd- This 
arises because Lboi spans a narrower range of values than 
Mbh- 

3. At a fixed value of Mbh, Lboi/LEdd increases systemati- 
cally along the sequence absorption-line nuclei — > tran- 
sition objects — > LINERs — > low-luminosity Seyferts — > 
high-luminosity Seyferts and quasars. There is consid- 
erable overlap among the classes. The apparent gap in 
Lboi/LEdd between the Palomar and SDSS sample may 
be an artifact of observational selection effects; the two 
surveys have very different sensitivity limits (Ho 2008). 

4. All emission-line nuclei in nearby galaxies are sub- 
Eddington systems, with the vast majority having 
Lboi/LEdd « 1. All LINERs and transition nuclei are 
characterized by Lboi/LEdd < 10" 2 . 

5. The combined distribution of L DO i or Lboi/Lfidd for the 
Palomar sample shows no evidence for bimodality or 
other indications of an abrupt transition between low- 
ionization (LINERs and transition objects) and high- 
ionization (Seyferts) sources. 

6. Because L^i spans a much larger range than M B h, Lboi 
broadly increases with increasing Lboi/LEdd- Low lu- 
minosity generally corresponds to low Eddington ra- 
tios. But there are important exceptions. A minor- 
ity of AGNs have low luminosities because they have 
low BH masses, not necessarily low Eddington ratios. 
Within the Palomar sample, NGC 4395 (Filippenko & 
Ho 2003) provides a good example (Fig. 4b), and sim- 
ilar types of low-mass AGNs have been discovered in 
SDSS (Greene & Ho 2004, 2007b). 
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FIG. 4. — Distribution of (a) BH masses and (b) Lbol vs. Z-bol/^Edd for objects separated by spectral classification. The bolometric luminosity is based on the 
2-10 keV X-ray luminosity, assuming L D0 i = 15.8Lx- The symbols are identified in the legend. The objects marked as "QSO" refer to the sample of high-luminosity 
Seyfert 1 nuclei and quasars studied by Greene & Ho (2007a). Line segments denote upper limits. 



4. SOURCES OF FUEL 

In this section we make some rough estimates of the mini- 
mum amount of fuel likely to be available in the central regions 
of nearby galaxies. For the moment, let us neglect any contribu- 
tion due to dissipation from a large-scale disk, external acqui- 
sition from tidal interactions and infall, or to discrete, episodic 
events such as tidal disruptions of stars. Galactic nuclei can be 
fed, in a steady state manner from the inner bulge of the galaxy, 
by at least two sources: (1) ordinary mass loss from evolved 
stars and (2) gravitational capture of gas from the hot interstel- 
lar medium. 

4. 1 . Stellar Mass Loss 

Present-day elliptical galaxies and the bulges of SOs and spi- 
rals contain mostly old, evolved stars. Red giants and planetary 
nebulae return a significant fraction of their mass to the inter- 
stellar medium through mass loss. For a Salpeter stellar initial 
mass function with a lower mass cutoff of 0.1 Mq, an upper 
mass cutoff of 100 M@, solar metallicities, and an age of 15 
Gyr, Padovani & Matteucci (1993) estimate 

M* »> 3 x 10" 11 ( ] Moyr" 1 . (2) 

This result is consistent, within a factor of ~2, with the work 
of Faber & Gallagher (1976), Ciotti et al. (1991), and Jung- 
wiert et al. (2001). Athey et al. (2002) obtained mid-infrared 
observations to probe more directly the mass-losing stars in el- 
liptical galaxies. They findM* w7x 10" 12 (L/L 0! b) M© yr" 1 . 
For B-V w 1 mag, appropriate for an evolved population (e.g., 
Fukugitaet al. 1995),M» w 2 x 10" 11 (L/L Qy ) M© yr" 1 , again 
close to the estimate by Padovani & Matteucci (1993). Thus, we 
use Padovani & Matteucci's relation to convert between V-band 
luminosity and mass loss rate. 



Hubble Space Telescope (HST) images have offered an un- 
precedently detailed view of the morphological structure of the 
central regions of nearby galaxies. The majority of galaxies 
contain central density concentrations, either in the form of nu- 
clear cusps or photometrically distinct, compact stellar nuclei 
(e.g., Lauer et al. 1995; Phillips et al. 1996; Carollo et al. 1997; 
Faber et al. 1997; Rest et al. 2001; Ravindranath et al. 2001; 
Boker et al. 2002; Ferrarese et al. 2006; Kormendy et al. 2009). 
The cusp profiles continue to rise to the resolution limit of HST 
(0"1), which is ~10 pc at a distance of 20 Mpc. The nuclear 
stellar population in most instances is old (Ho et al. 2003a; 
Sarzi et al. 2005; Zhang et al. 2008). How much gaseous ma- 
terial is available through stellar mass loss? We note that the 
uncertainties associated with the effectiveness of angular mo- 
mentum transport on large (1-10 kpc) or intermediate (0.1-1 
kpc) scales are bypassed by focusing only on nuclear (<10 pc) 
scales. Although the fate of the nuclear gas is not entirely clear, 
it is important to recognize that stellar mass loss confined to 
the nuclear cusp or nuclear cluster does furnish a steady state, 
in situ supply of gas that is in principle available for accretion. 
Shortly after being shed, the stellar gaseous envelopes quickly 
become thermalized with the hot ambient medium of the bulge, 
but some of the gas remains cool (Parriott & Bregman 2008). 
Even after the stellar ejecta joins the hot phase, the cooling time 
is sufficiently short in the inner region of the bulge that a cool- 
ing flow should develop (Mathews 1990). 

The three main Local Group galaxies (M31, M32, and M33) 
serve as instructive examples. From the work of Lauer et 
al. (1998), the central stellar densities of all three galax- 
ies rise steeply toward the center as p oc r" 1 ' 5±a5 ; at r = 0.1 
pc, the density reaches p w io 65±03 M© pc" 3 . More typi- 
cally, for galaxies beyond the Local Group, HST data probe 
r « 10 pc, where p pa 10— 10 3 L@y pc~ 3 for the "core" el- 
lipticals and p w 10 2 - 10 4 L©.y pc~ 3 for the "power-law" el- 
lipticals and bulges of early-type spirals and SOs (e.g., Faber 
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et al. 1997). Within a spherical region of r = 10 pc, the 
diffuse cores have Lw 4 x 10 4 -4 x 10 6 L Q y, which yields 
M* ra 1 x 10" 6 - 1 x 10" 4 M yr" 1 ; for the denser power-law 
cusps, Lw4 x 10 5 - 4 x 10 7 L Q . V , or ra 1 x 10" 5 -1 x 10" 3 
Mq yr" 1 . Centrally dominant nuclear star clusters, present in 
a large fraction of disk galaxies, typically have luminosities 
L ra 10 7 L (Lauer et al. 1995; Carollo et al. 1997; Boker 
et al. 2002), and hence M* ra 10" 3 M yr" 1 . 

4.2. Bondi Accretion 

The inner regions of ellipticals and bulges contain X-ray- 
emitting plasma, sustained through thermalized ejecta from 
stellar mass loss (Mathews 1990), with temperatures character- 
istic of the virial velocities of the stars, ~ 10 6 - 10 7 K. This dif- 
fuse, hot gas holds another plentiful fuel reservoir for accretion, 
through the mechanism described by Bondi (1952). The Bondi 
accretion rate depends on the gas density and temperature at 
the accretion radius, R a ra GM/c 2 , where c s ra O.U 1/2 km s" 1 
is the sound speed of the gas. From the continuity equation, 
Mb = 4-irR 2 a p a c s , where p a is the gas density at R a . Expressed 
in terms of typical observed parameters (see below), 



Mr 



7.3 x 10 



-4 



( M B 



\ 1O 8 M 



0.1 cm" 3 7 \ 



/ 200km s" 



M 



©yr 



(3) 

Recent high-resolution X-ray observations using Chandra 
find that the diffuse gas in the central few hundred parsec re- 
gions of elliptical galaxies typically has temperatures of kT ra 
0.3-1 keV and densities of n ra 0.1-0.3 cm" 3 (Di Matteo et 
al. 2001; Loewenstein et al. 2001; Sarazin et al. 2000; Pelle- 
grini 2005; Soria et al. 2006). Data for the bulges of spiral and 
SO galaxies are more fragmentary. Chandra has so far resolved 
the hot gas in the centers of a handful of bulges, including the 
Milky Way (Sbc; Baganoff et al. 2003), M81 (Sab; Swartz et 
al. 2003), NGC 1291 (Sa; Irwin et al. 2002), and NGC 1553 
(SO; Blanton et al. 2001). The center of M31 (Sb) has been 
investigated using both XMM-Newton (Shirey et al. 2001) and 
Chandra (Garcia et al. 2005). These studies suggest that bulges 
typically have gas temperatures of kT ra 0.3-0.5 keV. Informa- 
tion on gas densities is sketchier, but judging from the work on 
M81 and NGC 1291, a fiducial value might be n ra 0.1 cm" 3 . 

If, for simplicity, we assume that the hot gas in most bulges 
is characterized by n = 0.1 cm" 3 and kT = 0.3 keV, then 
M B ra 10" 5 - 10" 3 M Q yr" 1 for M BH = 10 7 - 10 s M Q . In ellip- 
tical galaxies M B h ra 10 8 - 10 9 M Q , and for n = 0.2 cm" 3 and 
kT = 0.5 keV, M B ra 10" 4 - 10" 2 M yr" 1 . We note that these 
estimates of the Bondi accretion rates are probably lower lim- 
its because the actual densities near R a are likely to be higher 
than we assumed. For the above fiducial temperatures and BH 
masses, R a ra 1 - 10 pc for bulges and ~ 10- 100 pc for ellip- 
ticals, roughly an order of magnitude smaller than the typical 
linear resolution achieved by Chandra for nearby galaxies. In 
the case of the Galactic center, n = 26 cm" 3 and kT = 1.3 keV 
at 10" (0.4 pc), rising to n = 130 cm" 3 and kT = 2 keV at 1" 
(Baganoff et al. 2003). 

4.3. Other Sources 



The two processes discussed above — ordinary stellar mass 
loss and Bondi accretion of hot gas — can be regarded as a 
conservative, steady state supply of fuel for galactic nuclei. 
Other sources, however, can raise this minimum level. In terms 
of purely stellar sources, some possibilities include (1) stellar 
mass loss enhanced by dynamical heating (Allen & Hughes 
1987; Armitage et al. 1996) or AGN irradiation (Edwards 1980; 
Scoville & Norman 1988; Voit & Shull 1988), (2) stellar-stellar 
collisions in a dense nuclear cluster (Spitzer & Saslaw 1966; 
Frank 1978; Rauch 1999), and (3) tidal disruption of stars by 
the central BH (Hills 1975; Rees 1988; Milosavljevic et al. 
2006). It is difficult to evaluate quantitatively the contribution 
these effects would make to the total fuel budget of nearby nu- 
clei; we merely note that cumulatively they may significantly 
boost the "baseline" accretion rate estimated above. 

We have also neglected any contribution from the cold phase 
of the interstellar medium. Nonaxisymmetric perturbations due 
to galaxy-galaxy tidal interactions, large-scale bars, nuclear 
bars, or nuclear spirals are often invoked as mechanisms for 
angular momentum transport of the cold gas in disk galaxies 
(e.g., Wada 2004). The effectiveness of these processes for fuel- 
ing nearby, relatively low-luminosity AGNs, however, has been 
unclear. With respect to the well-studied Palomar survey, AGN 
activity seems to be affected neither by large-scale bars (Ho et 
al. 1997c) nor by local galaxy environment (Schmitt 2001; Ho 
et al. 2003a). In any event, if dissipation of the cold gas does 
occur on nuclear scales, as inevitably it must at some level in at 
least some objects, it would further add to the fuel supply. 

5. IMPLICATIONS 
5.1. Accretion Flow and Radiative Efficiency 

The results presented in this paper provide some important 
insights into the nature of BH accretion in nearby galactic nu- 
clei. We have established, for the first time using a large, sta- 
tistically robust sample, that virtually all massive BHs in the 
nearby Universe share two common properties: they have low 
luminosities and radiate well below the Eddington limit. This 
holds for galaxies spanning a wide range of Hubble types and 
nuclear spectral classes. The median value of the bolometric 
luminosities are only L bo \ ra 10 39 - 10 41 erg s , and the median 
Eddington ratios range from Lboi/^Edd ra 2 x 10~ 7 to 3 x 10~ 4 . 

The extreme dimness of these nuclei strongly suggests that 
their accretion flows are radiatively inefficient. In the con- 
text of the class of accretion models commonly called opti- 
cally thin RIAF^] the accretion luminosity is given by (Ma- 
hadevan 1997) L acc = (77/O.I) [0.20(m/a 2 )] Me 2 , valid in the 
regime m > 10~ 3 a 2 , where M = mM^d and Msdd = 2.2 x 
10" 8 (77/O.I) (M B h/M ) M yr" 1 . This expression adopts the 
canonical values of the microphysics parameters used by Ma- 
hadevan (1997). In the notation used in this paper, 



0.7 



0.3 



/ Lbol \ 
V^Edd/ 



1/2 



(4) 



Thus, for L bol /L Edd = 10" 6 - 10" 4 , m ra 3 x 10" 4 -2 x 10" 2 , 
which lie comfortably within the regime of optically thin RI- 
AFs, m < m c rit ra a 2 ra 0.1 (Narayan et al. 1998). The corre- 
sponding absolute mass accretion rates, for Mbh ra 10 6 - 10 9 
M Q , are M ra 10" 5 - 10" 1 M yr" 1 . 

'As originally formulated by Narayan & Yi (1994, 1995), this class of accretion disk models was called advection-dominated accretion flows. Subsequent work 
has shown that such flows are inherently unstable to outflows and convection. To avoid delving into the technical details, which are unimportant for the present level 
of discussion, we simply follow Quataert (2001) and refer to this class of models as RIAFs. 
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The accretion rates estimated in § 4 provide a more direct 
argument that the radiative efficiency of the accretion flow, 
whatever its form, is likely to be low. The median bolomet- 
ric luminosities of the Palomar emission-line objects range 
from <~lxl0 40 erg s" 1 for transition objects to ~2xl0 41 erg 
s" 1 for Seyferts (Table 2). If this emission is produced by a 
canonical optically thick, physically thin disk, which radiates 
at L acc = rjMc 2 = 5.7 x 10 45 (r?/0.l) (M/M Q yr _1 ) erg s"\ we 
expect typical mass accretion rates of M « 2 x 10~ 6 to 4x 10~ 5 
M Q yr" 1 . These values of M are miniscule by comparison with 
the minimum accretion rates likely to be available through stel- 
lar mass loss and Bondi accretion alone. The majority of the 
Palomar AGNs are hosted by early-type disk galaxies (SOs and 
Sa-Sbc spirals; see Ho et al. 1997b, 2003a), whose bulges tend 
to have cuspy light profiles of the "power-law" type. As dis- 
cussed in § 4.1, the inner regions of such bulges should have 
M* « 10~ 5 - 10~ 3 Mq yr" 1 , and probably closer to the upper 
end of this range because of the additional contribution from 
central star clusters, which add M* « 10" 3 M yr" 1 . We have 
also erred on the side of caution by assuming that all of the 
stars are old. In actuality, the central regions of many spirals in 
the Palomar sample often show evidence for some contribution 
from composite populations (Ho et al. 2009), which will help 
to boost the mass loss rates even further. Bondi accretion of hot 
gas contributes roughly the same amount as stellar mass loss, 
M B ~ 10" 5 -10" 3 M Q yr" 1 . Thus, M tot = M* +M B « 10" 5 -10" 3 
Mq yr" 1 , but more likely, M tot ;> 10" 3 M Q yr" 1 . 

A similar exercise leads to an even stronger result for the 
absorption-line nuclei, which are found predominantly in el- 
liptical and SO galaxies. Here, the median Lboi of 2xl0 39 erg 
s" 1 requires only M«4x 10" 7 M Q yr" 1 for r\ = 0.1. On the 
other hand, the centers of the host galaxies can supply at least 
M* « 10" 6 - 10" 4 M Q yr" 1 for low-density cores, a factor of 
10 higher still in M* for power-law cusps, and yet another 10- 
fold increase for Mb. Ho et al. (2003b) highlighted the acute- 
ness of the luminosity-deficit problem for the nearest of the 
absorption-line objects, M32, whose 2.5 xlO 6 M Q BH emits 
merely 9.4x 10 35 erg s" 1 in the 2-10 keV band at an Eddington 
ratio of -3 x 10" 9 . 

These simple comparisons suggest that if r\ indeed is 0.1, 
then nearby galactic nuclei are 1-4 orders of magnitude under- 
luminous. For the accretion rates that we infer to be present, 
they should radiate far more prodigiously than actually ob- 
served. There are four possible interpretations of this finding. 
(1) First, our estimates of M tot could be too high by a large 
factor, namely 1-4 orders of magnitude. We consider this to 
be unlikely. Recall that M tot includes only normal mass loss 
from evolved stars and Bondi accretion of hot gas, either one 
of which alone would violate the luminosity limit for rj = 0. 1 ; 
we have conservatively neglected other potential sources of fuel 
(§ 4.3). (2) Second, it could be argued that perhaps the gas re- 
leased through stellar mass loss manages to escape from the nu- 
cleus before it gets accreted. In actively star-forming galaxies, 
for example, the collective effects of strong winds and shocks 
from massive stars can expel gas to large galactocentric dis- 
tances. This mechanism of gas removal, however, appears to be 
extremely implausible given the typical ages of the nuclear stel- 
lar population (Ho et al. 2003a). Storing the gas in an inert cold 
disk or converting it to young stars violates other observational 
constraints (Ho 2008). (3) Third, rj may be <C 0.1, as expected 
from RIAFs (Narayan et al. 1998; Quataert 2001). This is an 
argument that is frequently invoked to explain the apparent con- 



flict between the nuclear luminosities and Bondi accretion rates 
in some early-type galaxies (e.g., Fabian & Rees 1995; Ma- 
hadevan 1997; Di Matteo et al. 2000, 2001; Loewenstein et al. 
2001; Ho et al. 2003b). (4) Lastly, an inherent ambiguity exists 
between inefficient radiation and inefficient accretion. A va- 
riety of physical effects, summarized in Ho (2008), can divert 
the inflowing gas and severely curtail the amount of material 
that ultimately gets accreted. For instance, if RIAFs naturally 
develop outflows or winds, as recently shown in a number of 
studies (see Quataert 2001), the actual accretion rate would be 
much lower than that estimated at large radii. If so, then r\ may 
not need to be so exceptionally low, although it should still be 
substantially below 0. 1 because the outflow models ultimately 
rely on the accreting gas to be radiatively inefficient. Energy 
feedback from the AGN, associated with either disk outflows 
or small-scale radio jets — a ubiquitous feature of low-Lboi/^Edd 
systems (Ho 2002a) — may be another culprit for interrupting 
smooth mass inflow. Ho (2009) presents quantitative evidence 
that AGN feedback injects nongravitational perturbations to the 
kinematics of the ionized gas in the Palomar sources. 

5.2. Accretion States of Massive Black Holes 

By analogy with stellar BHs in X-ray binaries, supermas- 
sive BHs in galactic nuclei may evolve through different evo- 
lutionary phases, corresponding to distinct "states" (Narayan et 
al. 1998). The basic physical picture is that the structure of 
the accretion flow changes in response to variations in the ac- 
cretion rate. Parameterizing the accretion rate in terms of the 
dimensionless variable rh, one can define three, possibly four 
regimes. (1) When rh ;> 1, an object is in the "very high" 
state. The high radiation density in these "super-Eddington" 
sources traps the photons, and the accretion flow is that of an 
optically thick RIAF or slim disk (Begelman & Meier 1982; 
Abramowicz et al. 1988). An extragalactic analog of such sys- 
tems are the narrow-line Seyfert 1 nuclei (Pounds & Vaughan 
2000). (2) Objects satisfying m cr j t < rh < 1 correspond to those 
in the "high" state. These contain traditionally studied opti- 
cally thick, geometrically thin, radiatively efficient Shakura & 
Sunyaev (1973) disks, which are present in classical, relatively 
luminous Seyfert nuclei and quasars. (3) When rh < rh a \ u an 
optically thin, geometrically thick, radiatively inefficient flow 
develops, and the luminosity of the source plummets. Depend- 
ing on how low the accretion rate drops, one might distinguish 
between objects in the "low" state (10" 6 < rh < m CT n) versus 
those in the "quiescent" state {rh < 10 6 ). Objects in quiescence 
are those dominated by, or which exclusively contain, a pure 
RIAF, such as the Galactic center source Sgr A*. Low-state 
objects contain a hybrid structure consisting of an inner RIAF 
plus an outer thin disk, whose truncation radius recedes as rh 
decreases. Such a configuration has been suggested for a num- 
ber of low-luminosity AGNs, especially LINERs (Lasota et al. 
1996; Quataert et al. 1999; Ho et al. 2000; Ho 2002b, 2008). 
We note that the boundary between the low and quiescent states 
(m w 10~ 6 ) is purely illustrative; it remains to be demonstrated 
that there are two distinct states, and if so, where the transition 
truly lies. 

If AGN activity is characterized by distinct states, as 
schematically sketched above, then this ought to be reflected in 
the observed distribution of accretion rates for AGNs spanning 
the full range of rh. Since it is difficult to measure rh directly, 
Lboi/^Edd can be used as a surrogate for rh. We stress, however, 
that analysis of this kind is only meaningful when performed on 
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large, well-defined, statistically complete samples. It is danger- 
ous to combine samples with different selection criteria (e.g., 
Marchesini et al. 2004; Hopkins et al. 2006). With this in 
mind, we recall that the distribution of L\, i/L Edd for the entire 
Palomar sample (Figs, lb and 3b) — a more or less complete 
census of nearby galaxies — shows no obvious substructure that 
might be identified with physically distinct populations. Unfor- 
tunately, the volume sampled by the Palomar survey does not 
contain sufficient luminous AGNs to properly cover the upper 
end of the Lb i/L Edd distribution. This would be an important 
goal for future statistical studies of AGNs. 

6. SUMMARY 

Nearly all of the objects in the Palomar survey of nearby 
galaxies now have central stellar velocity dispersions, from 
which BH masses can be inferred using the Mbh~c* relation. 
A previously published collection of Ha line fluxes, in combi- 
nation with a newly assembled database of nuclear X-ray mea- 
surements and a reevaluation of the appropriate bolometric cor- 
rections, provides estimates of the accretion luminosity of the 
nuclei. We use these resources to evaluate the distribution of 
bolometric luminosities and Eddington ratios for a large, well- 
defined sample of galactic nuclei in order to investigate the na- 
ture of accretion onto massive BHs in nearby galaxies. 

Nearby galactic nuclei span at least 7 orders of magnitude in 
nuclear bolometric luminosities, fromLboi < 10 37 to ~ 3 x 10 44 
erg s" 1 , and an even broader range in Eddington ratios, from 
L bo i/L EAd « 10" 9 to 10" 1 . Both L b oi and L bol /L Edd , but espe- 
cially the latter, decrease systematically along the following 
spectral sequence: Seyferts — > LINERs — > transition objects 
— > absorption-line nuclei. The spectral diversity of emission- 
line nuclei reflects and is primarily controlled by variations in 
the mass accretion rate. The characteristic value of Lboi'^Edd 
also varies systematically along the Hubble sequence, increas- 
ing from galaxies with large to small bulge-to-disk ratios. 



The accretion rates inferred from the nuclear luminosities, 
assuming a standard radiative efficiency of rj = 0.1, are very 
low, typically M <; 10" 6 to 10" 5 M Q yr" 1 . Such tiny rates can 
easily be supplied in situ by ordinary mass loss from evolved 
stars in the nuclear stellar cusp or central star clusters, or by 
Bondi accretion of hot gas in the inner bulges of galaxies de- 
tected in X-ray observations. Indeed, we argue that conserva- 
tive estimates of the gas mass potentially available for accretion 
already far exceed the observational limits imposed by the lu- 
minosity measurements, suggesting that in many, if not most, 
nearby galaxies the radiative efficiency is likely to be much less 
than 0. 1 . The requirement for exceptionally low radiative ef- 
ficiencies, however, could be mitigated if the actual accretion 
rates are curtailed by AGN feedback in the form of disk out- 
flows or small-scale radio jets. The prevalence of RIAFs is fur- 
ther supported by the low Eddington ratios: all the objects in 
our sample are sub-Eddington systems, with the majority hav- 
ing values of Lboi'^Edd that satisfy the theoretically predicted 
criterion for RIAFs. 

We suggest that massive BHs in galactic nuclei evolve 
through distinct states in response to changes in the mass accre- 
tion rate. The nearby objects considered in this study are largely 
systems in the low or quiescent state. We see no evidence of bi- 
modality in the distribution of Lboi/^Edd, but this is probably a 
consequence of the limited volume probed by the Palomar sur- 
vey. It would be of considerable interest to extend the analysis 
presented here to include objects of higher luminosity in order 
to map out the full distribution of AGN luminosities and Ed- 
dington ratios. 
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Washington and by Chandra grant GO5-6107X. I thank Louis- 
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APPENDIX 

NUCLEAR X-RAY LUMINOSITIES 

This Appendix summarizes the X-ray luminosities used in this paper. We performed a comprehensive literature search of all 
published X-ray measurements for the AGNs (Seyferts, LINERs, and transition objects) and absorption-line nuclei in the Palomar 
survey. Because of the faintness of most nearby nuclei and potential confusion with circumnuclear emission, the most crucial 
consideration is angular resolution. With a point-spread function (PSF) of FWHM «1" and low background noise, the instrument 
of choice is ACIS on Chandra. Although less ideal, the HRI imager on ROSAT and the MOS camera on XMM-Newton, both having 
PSFs with FWHM 5", also yield acceptable data under most circumstances. For sources bright enough for rigorous spectral 
analysis, even lower resolution observations (e.g., ASCA) can be used if the nucleus can be isolated through spectral fitting. 

Of the 277 galaxies in the parent sample, acceptable literature data were located for 166. Most of the observations (75%) were 
acquired with Chandra/ACIS, and the rest were taken largely with ROSAT/HRl or XMM-Newton. Only a small handful come from 
ROSAT/PSPC and ASCA; one object was observed with BeppoSAX. We also performed a thorough search of the Chandra archives 
and included all useful, unpublished data that were nonproprietary as of 2007 October. A total of nine additional galaxies were 
located. We analyzed these data sets following standard techniques, as described in Ho et al. (2001) and Desroches & Ho (2009). 

Table 3 lists X-ray luminosities for the final sample of 175 galaxies. Because the literature data were acquired with a variety 
of different instruments and analyzed using many different techniques, we converted all the luminosities to one standard bandpass, 
2-10 keV. When reliable spectral fits are available, we use the published best-fit spectral slope to extrapolate to the desired bandpass. 
Otherwise, we assume a photon index of T = 1.8, which is close to the typical values observed in low-luminosity AGNs (see Ho 
2008, and references therein). Likewise, we quote intrinsic (absorption-corrected) luminosities whenever possible. Although many 
of the fainter sources do not have sufficient counts to constrain the absorbing column, many lines of evidence suggest that most 
low-luminosity AGNs do not suffer from much obscuration (Ho 2008). 
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TABLE 3 
Nuclear X-ray Luminosities 



Galaxy 


D L (Mpc) 


log (Lx/erg s x ) 


Tel. 


Ref. 


Galaxy 


D L (Mpc) 


log (Lx/erg s x ) 


Tel. 


Ref. 


(1)' 


(2) 


(3) 


(4) 


(5) 


(1)' 


(2) 


(3) 


(4) 


(5) 



IC 239 


16.8 


<38.05 


C 


1 


IC 1727 


8.2 


<37.92 


C 


2 


NGC 147 


0.7 


<36.18 


R 


3 


NGC 185 


0.7 


<35.92 


R 


4 


NGC 205 


0.7 


<36.41 


R 


3 


NGC 221 


0.7 


35.85 


C 


5 


NGC 224 


0.7 


35.85 


C 


6 


NGC 266 


62.4 


40.88 


C 


7 


NGC 315 


65.8 


41.63 


C 


8 


NGC 404 


2.4 


37.02 


c 


9 


NGC 410 


70.6 


40.09 


c 


1 


NGC 474 


32.5 


<38.36 


c 


1 


NGC 507 


65.7 


40.66 


c 


2 


NGC 628 


9.7 


38.14 


c 


2 


NGC 660 


11.8 


38.55 


c 


10 


NGC 676 


19.5 


39.00 


X 


11 


NGC 821 


23.2 


39.42 


c 


12 


NGC 1023 


10.5 


38.04 


c 


1 


NGC 1052 


17.8 


41.53 


c 


8 


NGC 1055 


12.6 


<38.19 


c 


13 


NGC 1058 


9.1 


<37.57 


c 


13 


NGC 1068 


14.4 


43.63 


A 


31 


NGC 1167 


65.3 


40.32 


c 


15 


NGC 1275 


70.1 


42.86 a 


X 


16 


NGC 1358 


53.6 


42.68 


B 


17 


NGC 1667 


61.2 


40.55 


c 


15 


NGC 1961 


53.1 


40.31 


R 


3 


NGC 2273 


28.4 


40.02 


R 


18 


NGC 2300 


31.0 


40.93 


R 


18 


NGC 2541 


10.6 


<38.00 


c 


19 


NGC 2639 


42.6 


40.85 


A 


20 


NGC 2655 


24.4 


41.87 


A 


20 


NGC 2681 


13.3 


38.72 


c 


8 


NGC 2683 


5.7 


38.01 


c 


19 


NGC 2685 


16.2 


39.90 


X 


11 


NGC 2768 


23.7 


40.59 


R 


21 



NGC 2775 


17.0 


39.07 


R 


18 


NGC 2787 


13.0 


38.79 


C 


13 


NGC 2832 


91.6 


<41.49 


R 


22 


NGC 2841 


12.0 


38.26 


C 


13 


NGC 3031 


1.4 


39.38 


C 


13 


NGC 3079 


20.4 


38.60 


c 


14 


NGC 3115 


6.7 


38.41 


c 


2 


NGC 3147 


40.9 


41.87 


c 


7 


NGC 3169 


19.7 


41.05 


c 


7 


NGC 3185 


21.3 


38.99 


X 


11 


NGC 3190 


22.4 


39.54 


R 


3 


NGC 3193 


23.2 


<39.41 


R 


3 


NGC 3226 


23.4 


39.99 


C 


8 


NGC 3227 


20.6 


41.70 


X 


11 


NGC 3245 


22.2 


39.29 


c 


8 


NGC 3368 


8.1 


39.10 


c 


14 


NGC 3377 


8.1 


38.24 


c 


23 


NGC 3379 


8.1 


37.53 


c 


8 


NGC 3384 


8.1 


37.79 


c 


3 


NGC 3412 


8.1 


<37.22 


c 


1 


NGC 3486 


7.4 


37.51 


c 


13 


NGC 3489 


6.4 


37.67 


c 


13 


NGC 3507 


19.8 


38.30 


c 


8 


NGC 3516 


38.9 


42.39 


X 


16 


NGC 3607 


19.9 


38.63 


c 


8 


NGC 3608 


23.4 


38.85 


c 


8 


NGC 3610 


29.2 


38.79 


c 


1 


NGC 3623 


7.3 


38.25 


c 


14 


NGC 3627 


6.6 


<37.22 


c 


13 


NGC 3628 


7.7 


38.24 


c 


8 


NGC 3640 


24.2 


<38.37 


c 


1 


NGC 3642 


27.5 


39.84 


R 


18 


NGC 3675 


12.8 


<37.98 


c 


13 


NGC 3718 


17.0 


40.44 


c 


19 


NGC 3884 


91.6 


41.89 


c 


24 


NGC 3898 


21.9 


<38.82 


c 


21 
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TABLE 3— Continued 



Galaxy 


D L (Mpc) 


log (Lx/erg s 1 ) 


Tel. 


Ref. 


Galaxy 


D L (Mpc) 


log (Lx/erg s 1 ) 


Tel. 


Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(1) 


(2) 


(3) 


(4) 


(5) 


NGC 3941 


18.9 


39.27 


X 


11 


NGC 4435 


16.8 


39.60 


R 


18 


NGC 3982 


17.0 


38.76 


c 


16 


NGC 4438 


16.8 


39.21 


C 


8 


NGC 3998 


21.6 


41.34 


R 


3 


NGC 4450 


16.8 


40.02 


C 


21 


NGC 4013 


17.0 


<37.57 


C 


10 


NGC 4457 


17.4 


39.05 


C 


8 


NGC 4036 


24.6 


39.96 


R 


18 


NGC 4459 


16.8 


38.87 


C 


8 


NGC 4051 


17.0 


42.07 


C 


14 


NGC 4472 


16.8 


<38.69 


C 


26 


NGC 4111 


17.0 


39.33 


C 


8 


NGC 4473 


16.8 


<38.20 


C 


1 


NGC 4125 


24.2 


38.94 


c 


8 


NGC 4477 


16.8 


39.60 


X 


11 


NGC 4138 


17.0 


40.11 


c 


11 


NGC 4478 


16.8 


<39.08 


c 


27 


NGC 4143 


17.0 


40.03 


c 


7 


NGC 4486 


16.8 


40.78" 


c 


8 


NGC 4150 


9.7 


<37.13 


c 


19 


NGC 4494 


9.7 


39.12 


c 


8 


NGC 4151 


20.3 


43.07 


c 


14 


NGC 4501 


16.8 


38.89 


c 


19 


NGC 4168 


16.8 


<38.24 


c 


16 


NGC 4527 


13.5 


38.80 


c 


10 


NGC 4192 


16.8 


39.58 


A 


20 


NGC 4548 


16.8 


39.74 


c 


7 


NGC 4203 


9.7 


39.69 


C 


13 


NGC 4550 


16.8 


<38.42 


c 


7 


NGC 4216 


16.8 


38.91 


c 


1 


NGC 4552 


16.8 


39.49 


c 


8 


NGC 4235 


35.1 


42.25 


X 


16 


NGC 4564 


16.8 


39.12 


c 


23 


NGC 4251 


9.7 


<37.78 


c 


1 


NGC 4565 


9.7 


39.56 


c 


7 


NGC 4258 


6.8 


40.89 


X 


34 


NGC 4569 


16.8 


39.41 


c 


13 


NGC 4261 


35.1 


40.59" 


c 


8 


NGC 4579 


16.8 


41.15 


c 


8 


NGC 4278 


9.7 


39.64 


c 


13 


NGC 4594 


20.0 


40.69 


c 


8 


NGC 4291 


29.4 


40.63 


R 


18 


NGC 4596 


16.8 


38.65 


c 


8 


NGC 4293 


17.0 


<39.37 


c 


4 


NGC 4621 


16.8 


38.85 


c 


28 


NGC 4314 


9.7 


38.13 


c 


8 


NGC 4636 


17.0 


39.38 


c 


8 


NGC 4321 


16.8 


<38.62 


c 


13 


NGC 4638 


16.8 


<39.11 


R 


3 


NGC 4346 


17.0 


<39.82 


c 


4 


NGC 4639 


16.8 


40.18 


c 


13 


NGC 4350 


16.8 


38.81 


c 


10 


NGC 4649 


16.8 


38.10 


c 


29 


NGC 4365 


16.8 


38.00 


c 


25 


NGC 4651 


16.8 


39.32 


R 


18 


NGC 4374 


16.8 


39.50" 


c 


8 


NGC 4698 


16.8 


38.69 


c 


8 


NGC 4382 


16.8 


37.87 


c 


25 


NGC 4713 


17.9 


38.40 


c 


2 


NGC 4388 


16.8 


42.14 


A 


33 


NGC 4725 


12.4 


39.11 


c 


13 


NGC 4394 


16.8 


<38.49 


R 


4 


NGC 4736 


4.3 


38.48 


c 


8 


NGC 4395 


3.6 


39.58 


C 


8 


NGC 4750 


26.1 


39.94 


c 


10 


NGC 4406 


16.8 


40.39 


R 


18 


NGC 4762 


16.8 


38.26 


c 


4 


NGC 4414 


9.7 


38.62 


C 


2 


NGC 4772 


16.3 


39.30 


c 


4 


NGC 4419 


16.8 


39.08 


c 


10 


NGC 4826 


4.1 


<37.33 


c 


13 



ACCRETION IN NEARBY GALAXIES 



TABLE 3— Continued 



Galaxy 


D L (Mpc) 


log (Lx/erg s 1 ) 


Tel. 


Ref. 


Galaxy 


D L (Mpc) 


log (Lx/erg s 1 ) 


Tel. 


Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(1) 


(2) 


(3) 


(4) 


(5) 


NGC 5005 


21.3 


39.94 


c 


10 


NGC 5850 


28.5 


<39.54 


R 


3 


NGC 5033 


18.7 


40.70 


c 


13 


NGC 5866 


15.3 


38.60 


C 


8 


NGC 5055 


7.2 


38.37 


c 


8 


NGC 6482 


52.3 


39.36 


C 


8 


NGC 5194 


7.7 


41.03 


B 


32 


NGC 6500 


39.7 


39.73 


C 


7 


NGC 5195 


9.3 


<38.00 


C 


13 


NGC 6503 


6.1 


37.10 


c 


14 


NGC 5273 


21.3 


41.61 


X 


11 


NGC 6654 


29.5 


<39.56 


R 


3 


NGC 5322 


31.6 


40.26 


R 


18 


NGC 7217 


16.0 


39.87 


R 


18 


NGC 5353 


37.8 


<39.29 


C 


18 


NGC 7331 


14.3 


38.66 


C 


8 


NGC 5354 


32.8 


<39.63 


R 


3 


NGC 7332 


18.2 


<39.49 


c 


27 


NGC 5363 


22.4 


<40.05 


X 


30 


NGC 7457 


12.3 


37.86 


c 


1 


NGC 5371 


37.8 


40.84 


R 


21 


NGC 7479 


32.4 


41.15 


X 


16 


NGC 5548 


67.0 


43.23 


X 


16 


NGC 7619 


50.7 


40.80 


c 


2 


NGC 5678 


35.6 


39.29 


C 


2 


NGC 7626 


45.6 


40.97 


c 


2 


NGC 5746 


29.4 


40.28 


C 


8 


NGC 7742 


22.2 


39.94 


R 


18 


NGC 5813 


28.5 


38.51 


c 


27 


NGC 7743 


24.4 


39.71 


X 


16 


NGC 5846 


28.5 


39.65 


c 


8 













NOTE. — Col. (1) Galaxy name. Col. (2) Adopted distance. Col. (3) X-ray luminosity in the 2-10 keV band. Col. (4) Telescope: 
A = ASCA; B = BeppoSAX; C = Chandra/ ACIS; R = ROSAT/HB1; X = XMM-Newton. Col. (9) Reference for the X-ray data: 
(1) Desroches & Ho (2009); (2) Chandra archives, analyzed in this paper; (3) Roberts & Warwick (2000); (4) Halderson et al. 
(2001); (5) Ho et al. (2003b); (6) Garcia et al. (2005); (7) Terashima & Wilson (2003); (8) Gonzalez-Martm et al. (2006); (9) 
Eracleous et al. (2002); (10) Dudik et al. (2005); (11) Cappi et al. (2006); (12) Pellegrini et al. (2007); (13) Ho et al. (2001); (14) 
Satyapal et al. (2004); (15) Pappa et al. (2001); (16) Panessa et al. (2006); (17) Deluit & Courvoisier (2003); (18) Liu & Bregman 
(2005); (19) Satyapal et al. (2005); (20) Terashima et al. (2002); (21) Komossa et al. (1999); (22) Dahlem fc Stuhrmann (1998); 
(23) Soria et al. (2006); (24) Pfefferkorn et al. (2001); (25) Sivakoff et al. (2003); (26) Loewenstein et al. (2001); (27) Pellegrini 
(2005); (28) Wrobel et al. (2008); (29) Soldatenkov et al. (2003); (30) Sansom et al. (2006); (31) Iwasawa et al. (1997a); (32) 
Fukazawa et al. (2001); (33) Iwasawa et al. (1997b); (34) Reynolds et al. (2009). 

"Nuclear luminosity possibly contaminated by jet emission. 



